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SUMMARY 


The  directional  damping  and  stability  of  a  fuselage  -  vertical- tail 
model  oscillating  freely  in  yaw  were  measured  at  a  Mach  number  of  O.l^ 
and  compared  with  the  damping  and  stability  obtained  by  consideration 
of  the  effects  of  unsteady  lift.  The  contribution  of  the  vertical  tail 
to  the  damping  in  yaw  of  the  model  agreed  within  the  limits  of  experi¬ 
mental  accuracy  with  predictions  of  the  approximate  finite-aspect-ratio 
unsteady-lift  theories  for  the  range  of  frequencies  tested,  and  the  two- 
dimensional  unsteady-lift  theory  predicted  a  much  greater  loss  in 
damping  with  reduction  in  frequency  than  was  shown  by  experiment  or  by 
the  approximate  finite-aspect-ratio  unsteady-lift  theories.  For  fre¬ 
quencies  comparable  to  those  of  lateral  airplane  motions,  both  the 
unsteady- lift  theories  and  the  experimental  results  indicated  that  the 
contribution  of  the  vertical  tail  to  the  directional  stability  of  the 
model  was  relatively  independent  of  frequency. 


INTRODUCTION 


The  advent  of  high-speed  airplanes  of  high  relative  density  has 
focused  attention  on  certain  problems  associated  with  the  dynamic  sta¬ 
bility  of  airplanes  which,  because  of  previous  unimportance,  have  here¬ 
tofore  been  neglected.  Among  these  problems  is  the  effect  of  the  perio¬ 
dicity  of  the  airplane  motion  on  the  effective  values  of  the  various 
stability  derivatives.  Reference  1  and,  more  recently,  references  2 
and  3  have  indicated  the  possibility  of  sizeable  effects  from  this 
problem. 

An  appreciable  amount  of  theoretical  work  on  these  unsteady-lift 
effects  exists  at  present;  however,  only  a  small  amount  of  experimental 
substantiation  of  the  results  is  available  for  the  low-frequency  range 
of  oscillation.  As  a  result  a  program  has  been  undertaken  in  the 
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Langley  stability  tunnel  to  determine  the  effects  of  such  variables  as 
frequency  and  amplitude  of  motion  on  the  contribution  of  the  various 
airplane  components  to  the  stability  derivatives  of  present-day  air¬ 
plane  configurations.  The  work  reported  herein  covers  that  phase  of 
the  investigation  which  considers  frequency  effects  on  the  directional 
damping  and  stability  of  a  model  undergoing  a  freely  damped  oscillatory 
yawing  motion.  The  effects  of  vertical-tail  aspect  ratio  and  compres¬ 
sibility  as  predicted  by  the  theoretical  treatments  are  discussed  in 
relation  to  the  experimental  stability  characteristics  obtained  by  the 
free-oscillation  and  by  the  curved-flow  procedures. 


SYMBOLS 


The  data  are  referred  to  the  system  of  stability  axes  and  are 
presented  in  the  form  of  standard  NACA  coefficients  of  forces  and 
moments  about  the  quarter-chord  point  of  the  mean  aerodynamic  chord  of 
the  normal  wing  location  of  the  model  tested.  (See  fig.  1.)  The 
coefficients  and  symbols  used  herein  are  defined  as  follows: 

Cl  lift  coefficient  (L/qSv) 

Cp  drag  coefficient  (D/qSy) 

Cy  lateral-force  coefficient  (Y/qSy) 

C2  rolling-moment  coefficient  (L'/qE^bw) 

Cm  pitching-moment  coefficient  (M/qSw‘Cw) 

Cn  yawing- moment  coefficient  (N/qSyby) 

L  lift,  pounds 

D  drag,  pounds 

Y  lateral  force,  pounds 

L'  rolling  moment,  foot-pounds 

M  pitching  moment,  foot-pounds 

N  yawing  moment,  foot-pounds 

q  dynamic  pressure,  pounds  per  square  foot 


NACA  TN  2657 


3 


¥ 


P 

V 

s 

b 

c 

ct 


c 


y 

A 


a 

V  P 


% 


t 

rb 

2V 

♦ 


r,i 


r,+ 


a) 

f 

k 


N 


* 


mass  density  of  air,  slugs  per  cubic  foot 

free-stream  velocity,  feet  per  second 

area,  square  feet 

span,  feet 

chord,  feet 

mean  tail  chord,  feet 


spanvise  distance  from  the  plane  of  symmetry 
aspect  ratio  (b2/s) 
angle  of  attack,  degrees 
angle  of  sideslip,  radians 


time,  seconds 

yawing- velocity  parameter 

angle  of  yaw,  radians 

yawing  velocity,  radians  per  second  (d\|r/dt) 
yawing  acceleration,  radians  per  second  per  second 

(a2*/at2) 

angular  velocity,  radians  per  second 
frequency,  cycles  per  second 

reduced-frequency  parameter  referred  to  semichord  of 
vertical  tail  (o>ct/2V) 

mechanical  spring  constant,  foot-pounds  per  radian 
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XZ 

tl/2 


m 


yawing  moment  of  inertia,  foot-pound.- second4 
time  to  damp  to  one-half  amplitude,  seconds 

damping  constant 

period  of  oscillation,  seconds 

nondimens ional  tail  length  referred  to  semichord  of 
vertical  tail 


(4) 


Z 

C(k)  =  F  +  iG 
F,G 

P  =  F  +  iG 
Q  =  H  +  iJ 
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kv 
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distance  from  origin  of  axes  to  midchord  point  of  vertical 
tail,  feet 


circulation  functions  used  by  Theodorsen  (reference  6) 

finite-span  functions  used  by  Biot  and  Boehnlein 
(reference  9) 

finite- span  correction  used  by  Reissner  (reference  8) 

compressible -flow  functions  used  in  place  of  F  and  G 
(references  10  and  11) 

functions  of  k  defined  in  text 
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Subscripts: 


wing 

vertical  tail 

amplitude 

due  to  friction 
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APPARATUS 

Model  and  Oscillating  Strut 


The  model  used  in  this  investigation  and  shown  in  figure  2  had 
^  separable  wing,  fuselage,  and  tail  surfaces  and  was  constructed  almost 
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entirely  of  balsa  wood  in  order  to  minimize  the  mass  and  make  the 
natural  frequency  of  the  model  on  its  mounting  as  high  as  possible. 

The  configuration  is  representative  of  that  of  several  present-day  high¬ 
speed  airplanes.  The  aspect  ratio  A  and  mean  chord  ct  of  the 
vertical  tail  were  2  and  4.9  inches,  respectively. 

The  model  was  mounted  at  the  quarter  chord  of  the  wing  mean  aero¬ 
dynamic  chord  on  a  1- inch-diameter  rod  which  rotated  with  respect  to  a 
circular  hollow  strut  and  which  was  connected  to  it  by  means  of  flexure 
pivots  made  of  Swedish  steel.  Flexure  pivots  were  used  rather  than 
bearings  in  order  to  minimize  the  friction  in  the  system.  The  pivots 
were  alined  with  the  center  of  the  rod  on  which  the  model  was  mounted. 
The  entire  assembly,  which  is  shown  in  figures  3  and  4,  was  fastened  to 
the  ceiling  of  the  tunnel  test  section  and  covered  by  a  streamlined 
fairing . 

The  inner  rod  or  oscillating  strut,  which  was  permitted  to  rotate 
and  upon  which  the  model  was  mounted,  extended  through  the  ceiling  of 
the  tunnel  test  section  and  was  capped  by  a  horizontal  crosspiece.  The 
period  of  oscillation  of  the  model  was  controlled  by  varying  the  moment 
of  inertia  about  the  vertical  axis  of  the  system.  This  variation  was 
accomplished  by  using  crosspieces  of  different  lengths  and  by  clamping 
weights  at  the  ends  of  the  crosspieces. 


Recording  Device 

A  continuous  record  of  the  displacement  in  yaw  of  the  model  was 
made  on  photographic  film  by  means  of  an  optical  recording  system.  A 
mirror  clamped  to  the  section  of  the  oscillating  strut  which  extended 
outside  the  tunnel  reflected  a  beam  of  light  from  a  point  light  source 
within  the  recorder.  The  beam  was  focused  on  a  roll  of  photographic 
film  by  means  of  an  8-inch-focal-length  double-concave  lens.  The  film 
roll  was  driven  at  constant  speed  by  gearing  it  to  an  electric  motor.  A 
10-cycle-per- second  neon  timing  light  incorporated  in  the  recorder 
exposed  timing  lines  on  the  film  in  order  that  time  as  well  as  model 
displacement  could  be  measured  accurately  on  the  film  record. 


TESTS 

Oscillation  Tests 


I 

The  free-damping  tests  were  conducted  with  the  model  in  an  inverted 
position,  in  order  to  minimize  strut  interference  effects,  and  at  zero 
angle  of  attack  with  respect  to  the  fuselage  reference  line.  The  tests 
were  conducted  on  the  combination  of  fuselage  and  tail  surface  and  on 
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the  fuselage  alone  "because  preliminary  tests  had  shown  that  the  presence 
of  the  wing  had  no  apparent  effect  on  the  damping  at  zero  angle  of 
attack.  The  range  of  moments  of  inertia  and  frequencies  covered  are 
shown  together  with  the  test  results  in  table  I.  Results  of  runs  made 
with  wind  on  and  off  are  given  in  order  to  illustrate  the  magnitude 
of  the  tares  applied. 

When  tests  were  made  of  the  fuselage  alone,  a  torque  rod  was 
employed  to  provide  a  stable  restoring  moment  to  the  fuselage  which 
was  equal  to  the  aerodynamic  moment  provided  by  the  vertical  tail  when 
attached.  The  torque  rod  was  clamped  between  the  extended  oscillating 
strut  and  a  building  structural  member  as  shown  in  figure  3*  The 
spring  constant  of  the  system  was  6.8  foot-pounds  per  radian  for  the 
flexure  pivots  alone  and  32.8  foot-pounds  per  radian  for  the  combina¬ 
tion  of  the  flexure  pivots  and  the  torque  rod. 


Force  Tests 

The  steady- state  stability  characteristics  of  the  model  were  deter¬ 
mined  through  an  angle-of-attack  range  from  0°  to  8°  for  the  fuselage 
alone,  the  fuselage  and  tail  surfaces,  and  the  complete  configuration 
by  means  of  wind-tunnel  force  tests.  The  model  was  mounted  on  a  single¬ 
strut  support  and  all  forces  and  moments  were  measured  by  a  six-component 
balance  system.  The  steady-state  sideslipping  derivatives  were  deter¬ 
mined  through  a  sideslip  range  of  ±5°,  and  the  yawing  derivatives,  at 
values  of  the  yawing- velocity  parameter  rb/2V  of  0,  -0.0291,  -O.0616, 
and  -0.08l0  by  standard  stability-tunnel  curved-flow  testing  procedures. 

All  tests  were  conducted  in  the  6-  by  6-foot  test  section  of  the 
Langley  stability  tunnel  at  a  dynamic  pressure  of  24.9  pounds  per  square 
foot,  a  Mach  number  of  about  0.l4,  and  a  Reynolds  number  of  442,000 
based  on  the  wing  mean  aerodynamic  chord.  Standard  jet-boundary  correc¬ 
tions  have  been  applied  to  the  angle-of-attack,  the  drag-coefficient, 
and  the  rolling-moment-coefficient  data  obtained  from  the  force  tests. 

In  addition,  the  steady- state  yawing  data  have  been  corrected  for  the 
effects  of  the  lateral  static-pressure  gradient  peculiar  to  the  curved- 
flow  testing  procedure.  Jet-boundary  corrections  were  not  applied  to 
the  oscillation  tests  because  of  thei-r  insignificance  for  the  test 
conditions  involved.  See  reference  4. 


REDUCTION  OF  OSCILLATION  TEST  DATA 


From  the  continuous  film  record  taken  of  the  motion  of  the  model 
after  an  initial  displacement,  the  amplitudes  of  the  successive  cycles 


8 


MCA  TN  2657 


were  measured  and  plotted  to  a  logarithmic  scale  against  time.  Inasmuch 
as  the  damping  is  of  a  logarithmic  nature,  the  resulting  plot  is  a 
straight  line  from  which  may  he  read  the  time  for  the  motion  to  damp  to 
one -half  amplitude  ^1/2)  • 


For  a  system  having  a  single  degree  of  freedom  in  yaw  such  that 
for  the  stahility-axis  system  the  angle  of  sideslip  3  is  the  negative 
of  the  angle  of  yaw  ty,  the  equation  of  motion  can  he  expressed  as 


(xz  -  kvK  -  (Nt  •  NP  +  N*f)*  -  (N*  -  np)*  =  0  (1) 

where  and  refer  to  the  mechanical  friction  of  the  system  and 

the  flexure-pivot  spring  constant,  respectively.  The  solution  of 
equation  (l)  can  he  written  as 

\|r  =  e"mt(A  sin  2itft  +  B  cos  2«ft)  (2) 


where 


2*f 


P*  -  NJ3 

'N*  -  Np  + 

pz  - 

2(IZ  -  »V) 

m 


»*  - 


The  damping  constant  may  he  expressed  in  terms  of  the  time  to  damp  to 
one-half  amplitude  as 


m 


-0.693 

tl/2 


The  directional  stability  and  damping  of  the  system  may  then  he 
expressed  in  nondimensional  form  in  a  manner  similar  to  that  of 
reference  5  as 


C 


nJ3 
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c 


n. 


-g.T72(iz  -  Nff)v 
qSwbv2 


wind  off 


(4) 


In  order  to  find  the  damping  due  to  friction  in  the  system,  the  damping 
of  the  model  is  determined  at  zero  wind  velocity  with  an  equivalent 
mass  in  place  of  the  vertical  tail  and  a  torque  rod  of  sufficient 
strength  to  produce  the  proper  frequency  inserted  in  the  system.  The 
term  is  interpreted  in  terms  of  a  factor  added  to  Cn^  because 

this  particular  group  of  free-damping  tests  afforded  no  convenient 
means  for  its  separation.  For  the  frequencies  employed  in  these  tests 

should  he  small  in  comparison  with  Cnp. 

The  damping  derivatives  were  obtained  from  data  similar  to  the 
sample  records  of  model  motion  and  plots  of  logarithm  of  amplitude 
against  time  given  as  figure  5*  A  small  residual  yawing  motion  of  the 
model  attributable  to  turbulence  in  the  tunnel  air  stream  necessitated 
determining  the  time  to  damp  to  one -half  amplitude  by  the  slope  of  the 
amplitude  envelope  at  the  largest  angles  of  yaw  in  order  to  minimize 
the  error  in  determining  the  damping  derivatives.  This  source  is 
believed  to  have  resulted  in  no  more  than  a  5“Percen_t  uncertainty  in 
the  results. 


the  factor  Cn^k' 


2/bw 

W 


AMLYTICAL  CONSIDERATIONS 


Various  unsteady-lift  theories  were  employed  to  calculate  the 
contributions  of  the  vertical  tail  to  the  directional  stability  and 
damping  of  the  model  for  comparison  with  the  experimental  results. 
These  theories  included  the  incompressible  cases  covered  by  Theodorsen 
for  two  dimensions  and  by  Jones  (reference  7) >  Reissner  (reference  8), 
and  Biot  and  Boehnlein  (reference  9)  with  various  degrees  of  approxi¬ 
mation  for  finite  span.  The  two-dimensional  work  including  the  sub¬ 
sonic  effects  of  compressibility  as  given  by  Possio  and  Frazer  and 
Skan  is  also  employed.  A  short  summary  indicating  procedures  for 
obtaining  the  necessary  derivatives  from  these  theories  is  included 
here  for  completeness. 


Two-Dimensional  Derivatives 

The  work  of  Theodorsen  on  the  derivation  of  the  expressions  for 
the  unsteady  lift  and  moment  of  a  two-dimensional  surface  undergoing 
sinusoidal  oscillations  in  an  incompressible  fluid  is  given  in 
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reference  6.  From  this  source,  the  moment  per  unit  span  about  the  air¬ 
plane  center  of  gravity  of  a  vertical  tail  mounted  on  an  airplane  which 
is  performing  sideslipping  oscillations  can  be  shown  to  be 


0, ac+3 

N(P)  =  -Jtqct2 


W 


(a  +  |)c(k)p 


(5) 


A  similar  expression  for  an  airplane  performing  yawing  oscillations  is 


N(t)  =  -§  qct3 


(6) 


An  appendix  to  this  paper  gives  a  discussion  of  the  conversion  of  these 
results  from  reference  6. 


The  Theodorsen  function  C(k)  is  a  complex  circulation  function  of 
the  reduced  frequency  and  in  terms  of  its  real  and  imaginary  parts  is 
written 


C(k)  =  F  +  iG 

The  functions  F  and  G  may  be  computed  from  expressions  given  in 
reference  6  in  terms  of  Bessel  functions  of  the  first  and  second  kinds 
of  orders  zero  and  one  and  of  argument  k.  In  stability  practice  equa¬ 
tions  (5)  and  (6)  are  usually  interpreted  in  terms  of  derivatives  as 
follows: 


qS^b.. 

NO)  =  CnpqSvbwP  +  Cn-  — 3 


N(*)  =  c  qSwbw2  *  +  C  qSwbw3  * 
mv'  ^  2V  ¥  +  ^  ^v2  ¥ 


(7) 


If  the  sinusoidal  motions 
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are  employed  for  a  substitution  into  equations  (5),  (6)  and  (?)>  the 
real  and  imaginary  terms  may  he  equated  and  the  derivatives  Cng>  Cnp' 

Cnr,  and  Cn£.  determined.  The  results  of  such  an  analysis  would  he  as 

follows: 


The  velocity  terms  contain  only  the  F  term  of  the  Theodorsen  function 
whereas  the  acceleration  terms  contain  the  G  term  divided  hy  k.  The 
corresponding  side-force  derivatives  are  derived  in  the  appendix.  For 
the  free -oscillation  tests  conducted  for  this  paper,  the  motion  is  such 
that  the  angle  of  sideslip  is  the  negative  of  the  angle  of  yaw.  The 
damping  derivative  measured  therefore  is  the  difference  between  the 
derivatives  and  Cn*  determined  from  equation  (8) .  This  fact 

was  illustrated  previously  with  respect  to  the  reduction  of  the  oscil¬ 
lation  test  data.  In.  a  like  manner,  the  directional  stability  parameter 
measured  can  he  shown  to  he  a  combination  of  the  Cnp  and  Cn^  terms, 

specifically 

%  *  cn.k2|s) 

Combination  of  the  results  of  equations  (8)  to  obtain  the  factors 
calculated  for  comparison  with  the  oscillation  tests  gives 
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B  -  -  (a  -  |)k  .  (a  +  i)20  ♦  (a2  - 
A  =  ^a2  +  ^)k2  +  +  -ij 2F  +  ^a2  -  |j2kG 


Finite-Span  Derivatives 

The  aerodynamic -span  effect  is  considered  in  reference  7  for  wings 
of  aspect  ratios  of  6  and  3  by  correcting  the  aerodynamic  inertia  and 
the  angle  of  attack  of  the  infinite-span  surface.  An  approximation 
employed  in  this  analysis  concerning  motions  of  long  wave  length  may 
make  the  results  subject  to  question  for  values  of  k  near  zero.  The 
reference  presents  expressions  for  calculating  the  finite-span  values 
of  the  Theodorsen  function  C(k).  The  use  of  these  F  and  G  circu¬ 
lation  functions  in  equations  (8)  and  (9)  results  in  values  of  the 
various  derivatives  considered  for  aspect  ratios  of  6  and  3. 

The  aerodynamic- span  effect  is  considered  for  wings  of  arbitrary 
aspect  ratio  in  reference  8.  The  three-dimensional  effect  of  the 
finite  span  may  be  obtained  by  adding  a  correction  terra  a  to  the 
basic  Theodorsen  two-dimensional  function.  The  finite- span  func¬ 
tion  C(k)  +  a  is  employed  in  a  rather  lengthy  and  tedious  computation 
in  order  to  determine  the  components  of  the  unsteady  lift  for  each  of 
five  points  along  the  span  of  the  oscillating  surface.  The  various 
stability  derivatives  may  then  be  determined  by  a  graphical  integration 
of  the  in-phase  and  out-of-phase  components.  An  approximation  made  in 
this  analysis  is  of  such  a  nature  as  to  make  the  method  primarily 
suited  to  the  higher  aspect  ratios.  A  practical  lower  limit  of  appli¬ 
cability  is  not  known. 

A  one-point  approximation  to  the  span  effect  for  the  calculation 
of  the  unsteady  lift  may  be  employed  rather  than  the  five-point  proce¬ 
dure.  Comparative  calculations  made  for  the  present  investigation 
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have  shown  the  one-point  approximation  to  yield  results  for  low  aspect 
ratios  which  are  essentially  the  same  as  those  calculated  by  the 
longer  procedure. 

Biot  and  Boehnlein  (reference  9)  considered  the  aerodynamic- span 
effect  on  unsteady  lift-  by  an  approach  that  yields  a  closed  form  for  a 
one-point  approximation  to  the  span  effect.  The  finite-span  circula¬ 
tion  functions  used  are 

P  =  F  +  iG 


t 


Q  =  H  +  iJ 

For  infinite  aspect  ratio 

P  =  Q  =  C(k) 

where  C(k)  is  the  two-dimensional  Theodorsen  function.  The  func¬ 
tions  P  and  Q  may  be  calculated  for  the  argument  k  by  methods 
given  in  reference  9,  and  equations  (8)  and  (9)  are  again  used  to 
calculate  the  various  stability  derivatives.  The  factors  A  and  B  _ 
which  appear  in  equations  (9)  must,  however,  be  written  in  terms  of  P 
and  Q  as  follows: 


A  =  (2a  -  l)akG  +  ^a  -  |)kj  +  {a2  +  ^)k2  +  2elF  +  H 


B  =  (2a  -  l)akF  + 


Effect  of  Compressibility 

Subsonic  compressibility  effects  on  the  unsteady  damping  in  yaw 
and  directional  stability  of  a  two-dimensional  oscillating  lifting 
surface  can  be  considered  through  the  use  of  reference  10  which  is  an 
extension  of  the  calculations  of  Possio  (reference  11) .  The  infinite - 
span  compressible-flow  values  for  the  stability  derivatives  can  be 
obtained  by  the  use  of  equations  (8)  and  (9)  for  Mach  numbers  0,  0.5> 
and  0.7.  In  the  notation  of  Frazer  and  Skan  (reference  10) ,  however, 

A  =  ^3  "  t(t  +  ^(Ml  +  z3  -  I  Zl)  +  Ml 

B  =  +  *)(m2  +  "  1  Z2)  +  M2 
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The  correspondence  between  the  notation  of  Theodorsen  and  that  of 
Frazer  and  Skan  is  shown  in  reference  12.  The  values  of  the  Z  and 
M  functions  are  tabulated  by  Frazer  and  Skan  for  Mach  numbers  0,  0.5, 
and  0,7  and  a  significant  range  of  reduced  frequencies.  Reissner  has 
recently  extended  his  work  on  oscillating  wings  of  finite  span  to 
include,  to  the  same  degree  of  approximation  as  for  the  incompressible 
case,  the  subsonic  effects  of  compressibility.  (See  reference  13.) 


RESULTS  AND  DISCUSSION 
Presentation  of  Results 


The  experimental  values  of  the  longitudinal,  directional,  and 
yawing  characteristics  of  the  complete  model,  the  model  with  wing 
removed  and  the  fuselage  alone  as  determined  by  force  tests  are  given 
in  figure  6.  The  experimental  values  of  the  damping  in  yaw  of  the 
fuselage  and  tail  assembly  (designated  FT)  and  the  fuselage  alone 
(designated  F)  obtained  by  the  freely  damped  oscillation  technique  for 
several  frequencies  of  oscillation  are  shown  in  figures  7  and  8.  The 
frequency  and  aspect-ratio  effects  on  the  damping  in  yaw  contributed  by 
the  vertical  tail  as  calculated  from  the  unsteady-lift  theories  dis¬ 
cussed  previously  are  also  shown  in  these  figures.  Figure  9  gives  a 
comparison  of  the  relative  importance  of  Cnr  and  Cn £  on  the  damping 

in  yaw  of  the  model.  Figures  10  to  12  are  similar  figures  summarizing 
frequency  and  aspect-ratio  effects  on  directional  stability.  Figures  13 
and  Ik  give  the  computed  effects  of  compressibility  on  the  vertical- 
tail  contribution  to  the  damping  in  yaw  and  directional  stability  of 
the  test  model. 


Discussion 

The  variation  of  the  experimentally  determined  damping  in  yaw  of 
the  test  model  with  frequency  is  in  reasonably  good  agreement  for  the 
range  covered  by  the  tests  with  the  various  approximate  theoretical 
treatments  which  consider  the  effect  of  finite  span.  (See  figs.  7  and  8.) 
The  contribution  of  the  fuselage  to  the  damping  in  yaw  is  small  for  the 
range  of  test  frequencies.  A  two-dimensional  approximation  to  the 
damping  as  is  frequently  employed  for  flutter  work  appears  to  be 
remarkably  good  for  k  >  0.1  and  for  aspect  ratios  of  6  and  greater 
(fig.  8).  A  large  effect  of  aspect  ratio,  however,  is  shown  for  aspect 
ratios  below  6  for  the  entire  frequency  range  and  for  all  aspect  ratios 
for  frequencies  less  than  k  =  0.1.  The  two-dimensional  result  indi¬ 
cates  a  loss  in  damping  at  low  frequencies  much  greater  than  the  finite- 
span  result  and,  in  fact,  shows  a  reversal  in  sign  of  the  damping  at 
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sufficiently  low  frequencies.  A  number  of  papers  dealing  with  unsteady 
lift  have  noted  this  characteristic  of  the  two-dimensional  aerodynamic 
forces.  (See  references  1  to  3  and  l4.) 


An  examination  of  the  values  of  the  derivatives  Cn_^  and  Cn^  as 

computed  from  reference  6  for  infinite  aspect  ratio  and  reference  7  for 
aspect  ratios  of  3  and  6  (fig.  9)  indicates  that  the  derivative  Cn£  is 

critically  dependent  on  reduced  frequency  for  the  low  range  of  reduced 
frequencies  and  infinite  aspect  ratio.  The  derivative  Cn^  is  only 

slightly  dependent  on  the  reduced  frequency  for  infinite- aspect  ratio 
and  even  less  so  for  the  lower  aspect  ratios.  The  dependence  of  Cnp 

on  reduced  frequency  is  responsible  for  the  decrease  in  damping  shown 
for  the  infinite-aspect-ratio  case  as  the  frequency  is  reduced  (figs.  7 
and  8) . 


The  aspect-ratio-2  results  shown  in  figures  7  and  10  were  obtained 
from  the  theory  of  Jones  (reference  7)  for  aspect  ratio  3  by  use  of  a 
correction  factor  which  was  the  ratio  of  the  lifting-line-theory  lift- 
curve  slopes  for  aspect  ratios  of  2  and  3*  This  correction  is  only 
approximate  and  was  made  in  order  to  place  the  results  on  a  comparable 
basis.  The  variation  of  Cn^  -  Cn^  with  reduced  frequency,  shown  in 

figure  7,  indicates  no  appreciable  effect  of  frequency  to  zero  frequency, 
In  fact,  the  result  of  Jones  for  aspect  ratio  3  reduces  to  a  finite 
value  at  k  =  0  given  by 


Cnr  -  cnp 


(l.2a2 


0.370a  +  0.515) 


which  is  not  the  case  for  the  two-dimensional  result  that  approaches 
infinity  as  k  approaches  zero.  This  result  is  not  expected  to  be 
very  accurate  close  to  k  =  0,  however,  because  of  the  nature  of  the 
assumptions  made  in  its  derivation  as  mentioned  previously. 

When  the  finite- span  effects  are  considered,  it  appears  that  these 
effects  are  especially  significant  with  regard  to  the  damping  in  yaw  of 
this  model  for  aspect  ratios  smaller  than  6.  (See  fig.  8.)  The  theory 
also  indicates  a  significant  decrease  in  the  damping  at  low  frequencies 
when  the  aspect  ratio  was  increased  from  20  to  00. 

The  experimental  values  of  damping  in  yaw  obtained  by  free  oscil¬ 
lation  agree  closely  with  the  values  of  Cn^  obtained  by  the  standard 

curved-flow  testing  procedure  used  in  the  Langley  stability  tunnel.  The 
curved- flow  technique  results  in  the  k  =  0  value  of  Cn  .  No  unsteady 
effects,  however,  are  obtained  by  this  technique. 
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The  unsteady-lift  theories  indicate  that,  at  frequencies  of  oscil¬ 
lation  from  k  =  0  to  k  =  0.1,  the  frequency  has  little  effect  on  the 
directional  stability  contributed  by  the  vertical  tail  but  that  high 
frequencies  of  oscillation  can  cause  the  directional  stability  con¬ 
tributed  by  the  vertical  tail  to  decrease  greatly  and  eventually  to 
change  sign  (figs.  10  and  ll)  .  Changes  in  aspect  ratio,  for  low  aspect 
ratios,  have  a  large  effect  on  the  directional  stability;  for  high 
aspect  ratios,  changes  in  aspect  ratio  are  less  significant.  The 
experimental  values  of  directional  stability  obtained  from  the  free- 
damping  tests  remained  constant  for  the  range  of  frequencies  investi¬ 
gated,  although  the  magnitudes  are  somewhat  smaller  than  those  predicted 
for  the  tail  contribution  by  the  calculation  methods  for  aspect  ratio  2. 
This  result  is  not  unexpected  because  all  the  finite-span  unsteady-lift 
theories  employed  maintain  some  resemblance  to  the  Prandtl  lifting-line 
theory  which,  of  course,  yields  values  of  the  lift- curve  slope  for  low 
aspect  ratios  which  are  unrealistically  high.  The  values  of  directional 
stability  obtained  by  both  the  free-oscillation  and  curved-flow  proce¬ 
dures  were  about  the  same  .magnitude. 


Figure  12  gives  a  comparison  of  the  contribution  of  a  vertical  tail 
to  the  stability  derivatives  Cn^  and  Cn^  as  computed  from  refer¬ 
ences  6  and  7-  These  derivatives  were  combined  in  the  proportion 


Cn  +  Cn.k 
n0 


2  (hi 

nr“  yct 


I 


for  comparison  with  the  results  of  the  oscillation  tests.  Although  Cn. 

2  r 


is  of  appreciable  magnitude  itself,  the  factor  C 


k2(M  i 

n.K  1 -  1 

r  \Ct  / 


is  very 


small  relative  to  Cn^  for  the  range  of  frequencies  around  k  =  0.01 

and  thus  should  not  greatly  affect  the  directional  stability  of  the 
usual  airplane. 


The  computed  two-dimensional  Mach  number  effects  on  the  damping  in 
yaw  and  directional  stability  of  the  tail  surface  of  the  test  model  are 
shown  in  figures  13  and  l4.  These  results  were  included  to  illustrate 
the  powerful  effect  of  compressibility  on  the  damping  contributed  by 
the  tail  surface.  A  complete  coverage  of  one-degree-of-freedom  flutter 
by  use  of  these  aerodynamic  forces  which  may  be  obtained  from  refer¬ 
ences  10  to  12  is  given  in  reference  14.  Compressibility  effects 
decrease  the  tail-surface  damping  at  low  frequencies  and  increase  the 
damping  at  higher  frequencies.  The  directional  stability  increases  at 
all  frequencies  with  increase  in  Mach  number  with  the  largest  increases 
occurring  in  the  high-frequency  range  of  oscillations.  Reference  15 
which  presents  the  two-dimensional  oscillating  aerodynamic  forces  at 


t 


4 
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sonic  speed  indicates  larger  deviations  from  the  results  for  incompress¬ 
ible  flow  than  for  the  subsonic  results  given  here. 


CONCLUSIONS 


The  following  conclusions  were  drawn  from  the  results  of  an  inves¬ 
tigation  made  at  a  Mach  number  of  0.l4  of  the  effects  of  frequency  and 
aspect  ratio  on  the  directional  damping  and  stability  of  a  fuselage  and 
vertical  tail  oscillating  freely  in  yaw: 

1.  The  contribution  of  the  vertical  tail  to  the  damping  in  yaw  of 
the  model  agreed  within  the  limits  of  experimental  accuracy  with  predic¬ 
tions  of  the  approximate  finite-aspect- ratio  unsteady-lift  theories  for 
the  range  of  frequencies  tested. 

2.  Two-dimensional  unsteady-lift  theory  predicted  a  much  greater 
loss  in  damping  with  reduction  in  frequency  than  was  shown  by  experi¬ 
ment  or  by  the  approximate  finite-aspect-ratio  unsteady-lift  theories. 

3.  Both  the  unsTteady-lift  theories  and  the  experimental  results 
indicated  that  the  contribution  of  the  vertical  tail  to  the  directional 
stability  of  the  model  was  relatively  independent  of  frequency  for 
frequencies  comparable  to  those  of  lateral  airplane  motions. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  November  26,  1951 
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APPENDIX 

m 

CONVERSION  OF  UNSTEADY  FORCES  AND  MOMENTS  FROM 
FLUTTER  TO  STABILITY  AXES 


The  conversion  of  the  contribution  of  the  vertical  tail  to  the 
yawing  moment  of  an  airplane  in  terms  of  unsteady-lift  functions  and 
the  Theodorsen  flutter  axes  to  the  stability  axes  which  are  fixed  in 
the  body  is  included  here  as  a  matter  of  general  interest.  Reference  l6 
gives  considerable  attention  to  such  axis  conversions.  The  following 
sketches  indicate  the  paths  and  attitudes  of  airplanes  performing 
sinusoidal  sideslipping  and  yawing  oscillations  within  the  limits  of 
the  stability-axis  system: 


W 

* 


Yawing,  t  =  ^Qe 
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The  Theodorsen  coordinates  h  and  a  and  the  stability-axis  coordi¬ 
nates  3  and  i  are  shown  with  respect  to  the  paths  followed  hy  the 
airplanes  in  sideslipping  and  yawing  flight. 

In  sideslipping  flight 


a  =  0 


and  in  yawing  flight 


A 


From  reference  d  the  contribution  of  the  vertical  tail  to  the  yawing  moment 
of  an  airplane  in  terms  of  the  flutter  coordinates  is 


N(a,h) 


pet2  " 

'/1  \Vct 

\  anct  -I 

~T~  1 

L 

X  \2  "  a/  2 

a  +  k  ^  +  a 

^  1 

0 J 

d 

PV 


*(a  +  |)c(k)  jva  +  h  +  -7^  -  a)  ij 


If  this  moment  is  for  a  sideslipping  motion  or  a  yawing  motion  in  the 
sense  given  by  stability  axes  fixed  in  the  body  and  the  following  terms 
are  substituted  for  the  sideslipping  case 


a  -  0 
a  =  0 
a  =  0 

h  =  -VP 

■  •  • 

h  =  -V|3 


ms 


¥ 
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and  for  the  yawing  case 


=  -Vi) f 
h  =  -Vt 


a  = 

•  • 

a  =  ^ 

•  •  M 

a  = 

the  following  results  Eire  obtained  for  the  sideslipping  case 

2  \  2 
pc.  /ajtc.V  \  pVc .  ji.  n 

N(p)  - - Pj  + - £-(*  +  |jC(k)(-V3) 

and  for  the  yawing  case 
2 


NO)  =- 


Pct 


f/1  \Vct  •  rtct  /I  p\~  a*ctV  r 

Ki  '  *)■ -r  * +  -r-S  *  “1*  +  -T-  *J 


¥  + 


pVc. 


JI (a  +  |)c(k)fjr(t  -  a)ij 


which  reduce  for  the  sideslipping  case  to 

nO)  =  -*qct2jjw^  +  ^  +  |)c(k)p] 

and  for  the  yawing  case  to 


N(i)  =-|  qct3B  -  C(k)(i  -  a®)]*  +  +  ^4 


Similar  expressions  may  he  obtained  for  the  side  force  in  pure  side¬ 
slipping  and  yawing  flight.  Reference  6  gives  the  side  force  of  the 
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vertical  tail  in  terms  of  the  flutter  axes  as 

! 

2  ,  % 

Pct  /  .  itc^aaX  r-  .  ct/l 

Y(a,h)  =-— —  (vita  +  nh  -  — —  I  -  jrpVctC(k)  Va  +  h  +  —h; 

which,  when  the  following  terms  are  substituted  for  the  sideslipping 
case 

h  =  -V3 

>  •  • 

h  =  -Vf3 

a  =  0 
a  =  0 
a  =  0 

and  for  the  yawing  case 

h  =  -V)|r 
h  =  -Vt 

a  =  >|f 

•  • 

a  = 

a  =  V 

becomes  for  pure  sideslipping  and  yawing  flight 

Y(P)  =  qct  |[ct(|)  +  *tC(k) p| 

and 

Y(+)  --  qct2  |jc(k)(i  -  a)|  -  ^  a  Yj 
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Substituting  the  motions,  in  the  manner  done  in  the  text,  for  the 
yawing  moment, 


P  =  IV 


ia)t 


t  e' 
o 


icst 


into  these  expressions  and  into  the  expressions  in  terms  of  stability 
derivatives  which  are 

.  •  ^l^w^w 

Y(P)  =  Cy^PqS*  +  CY.0  -w- 


Y(i)  =  ~ 2V  + 


gives  the  following  relationships  for  the  tail  contribution  to  the  side 
force  derivatives  when  the  real  and  imaginary  parts  of  these  equations 
are  equated: 
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Figure  1.-  System  of  stability  axes.  Arrows  indicate  positive  forces, 
moments,  and  angular  displacements.  Yaw  reference  generally  chosen 
to  coincide  with  initial  relative  wind. 
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a  Recorder 


b  Lens 
c  Film 
d  Timing  light 

e  Light  source 


Figure  3.-  Diagram  of  test  equipment. 


Figure  Test  model  mounted  on  free -oscillation  strut  in  Langley 

stability  tunnel. 
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Figure  11.-  The  effect  of  aspect  ratio  on  Cnp  as  calculated  by  the 
method  of  Biot  for  a  range  of  reduced  frequencies. 


(b)  Range  of  reduced  frequencies. 

Figure  13®-  The  effect  of  compressiblity  on  the  damping  in  yaw 


